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Abstract We deal with the attempts to measure the Lense-Thirring effect with the 
Satellite Laser Ranging (SLR) technique applied to the existing LAGEOS and LA- 
GEOS II terrestrial satellites and to the recently approved LARES spacecraft. Accord- 
ing to general relativity, a central spinning body of mass M and angular momentum S 
like the Earth generates a gravitomagnetic field which induces small secular precessions 
of the orbit of a test particle geodesically moving around it. Extracting this signature 
from the data is a demanding task because of many classical orbital perturbations 
having the same pattern as the gravitomagnetic one, like those due to the centrifugal 
oblateness of the Earth which represents a major source of systematic bias. The first 
issue addressed here is: are the so far published evaluations of the systematic uncer- 
tainty induced by the bad knowledge of the even zonal harmonic coefficients of the 
multipolar expansion of the Earth's geopotential reliable and realistic? Our answer is 
negative. Indeed, if the differences AJg among the even zonals estimated in different 
Earth's gravity field global solutions from the dedicated GRACE mission are assumed 
for the uncertainties 5 Ji instead of using their covariance sigmas a j e , it turns out that 
the systematic uncertainty Sfi in the Lense-Thirring test with the nodes Q of LAGEOS 
and LAGEOS II may be up to 3 to 4 times larger than in the evaluations so far pub- 
lished (5 — 10%) based on the use of the sigmas of one model at a time separately. 
The second issue consists of the possibility of using a different approach in extracting 
the relativistic signature of interest from the LAGEOS-type data. The third issue is 
the possibility of reaching a realistic total accuracy of 1% with LAGEOS, LAGEOS II 
and LARES, which should be launched in November 2009 with a VEGA rocket. While 
LAGEOS and LAGEOS II fly at altitudes of about 6000 km, LARES will be likely 
placed at an altitude of 1450 km. Thus, it will be sensitive to much more even zonals 
than LAGEOS and LAGEOS II. Their corrupting impact has been evaluated with the 
standard Kaula's approach up to degree I — 60 by using AJg and aj e ; it turns out that 
it may be as large as some tens percent. The different orbit of LARES may also have 
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some consequences on the non-gravitational orbital perturbations affecting it which 
might further degrade the obtainable accuracy. 

Keywords Experimental tests of gravitational theories • Satellite orbits • Harmonics 
of the gravity potential field 

PACS 04.80.Cc • 91.10.Sp ■ 91.10.Qm 



1 Introduction 

In the weak-field and slow motion approximation, the Einstein field equations of general 
relativity get linearized resembling to the Maxwellian equations of electromagntism. 
As a consequence, a gravitomagnetic field B s , induced by the off-diagonal components 
goi,i = 1,2,3 of the space-time metric te nsor related to t he mass-energy currents of 
the source of the gravitational field, arises |Mashhoonll2007r i. The gravitomagnetic field 



affects orbiting test particles, pr ecessing gyroscopes, moving cl o cks and atoms and prop- 
agating electromagnetic waves (|Ruggiero and T artaglia 2002; Schafcr 2004). Perhaps, 
the most famous gravitom agnetic effects are the precession of the axis of a gyroscope 



dPughl Il959l ; ISchiflj Il96d ) and the Lense-Thirrin£] precessions |Lense and Thirrind 
1918) of the orbit of a test particle, both occurring in the field of a central slowly 
rotating mass like, e.g., our planet. Direct, undisputable measurements of such funda- 
mental predictions of general relativity are not yet available. 

The measurement of the gyroscope precession in the Ea rth's gravitational field has 



been the goal of the dedicated space-basec0 GP-B mission |Everittlll973 ;l Everitt et al.l 
200 if ) launched in 2004 and carrying onboard four superconducting gyroscopes; its data 
analysis is still ongoing. The target accuracy was originally 1%, but it is still unclear if 
the GP-B team will succeed in reaching such a goal because of some unmodelled effects 
affecting the gyroscopes: 1) a time variation in the polhode motion of the gyroscopes 
and 2) very large classical misalignment torques on the gyroscopes. 

In this Chapter we will focus on the attempts to measure of the Le nse-Thirring 
effect in the gravitation al fiel d of the Ear th; for Mars and the Sun see (jloriol 12006a; 
Kroghll2007l ; Ilorioll2007ah and (|lorioll2008ah , respectively. Far from a localized rotating 



body with angular momentum S the gravitomagnetic field can be written as 

B s = --^[S-3(S-f)f], (1) 

where G is the Newtonian gravitational constant and c is the speed of light in vacuum. 
It acts on a test particle orbiting with a velocity v with the non-central acceleration 
(|Soffelll 19891 ) 

A hT = --vxB s (2) 

c 

which induces secular precessions of the longitude of the ascending node SI 

^ = c2a 3 ( l_ e2) 3/2' ( 3 ) 



1 According to an interesting historical analysis recently performed by Pfistcr ( 20071') . it 
would be more correct to speak about an Einstein- Thirring-Lense effect. 

2 See on the WEB http: / /einstein. Stanfor d. edu/| 
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Fig. 1 Keplerian orbit. The longitude of the ascending node (2 is counted from a reference 
X direction in the equator of the central body, assumed as reference plane {X, Y}, to the line 
of the nodes which is the intersection of the orbital plane with the equatorial plane of the 
central body. It has mass M and proper angular momentum S. The argument of pericentre 
w is an angle in the orbital plane counted from the line of the nodes to the location of the 
pericentre, here marked with II. The time-dependent position of the moving test particle of 
mass m is given by the true anomaly /, counted anticlockwise from the pericentre's position. 
The inclination between the orbital and the equatorial planes is i. Thus, Q,w,i can be viewed 
as the three (constant) Euler angles fixing the configuration of a rigid body, i.e. the orbit which 
in the unperturbed Keplerian case does change neither its shape nor its size, in the inertial 
{X, Y,Z} space. Courtesy by H.I.M. Lichtenegger, IWF, Graz. 




and the argument of pericentre ui 



6GS cos i 



c 2 a 3 (l 



2)3/2 ' 



(4) 



of the orbit of a test particle. In eq. (J3j and eq. (|4} a and e are the semimajor axis 
and the eccentricity, respectively, of the test particle's orbit and i is its inclination to 
the central body's equator. The semimajor axis a fixes the size of the ellipse, while its 
shape is determined by the eccentricity < e < 1; an orbit with e = is a circle. The 
angles Q and u establish the orientation of the orbit in the inertial space and in the 
orbital plane, respectively. £2, to and i can be viewed as the three Euler angles which 
determine the orientation of a rigid body with respect to an inertial frame. In Figure 
Q]we illustrate the geometry of a Keplerian orbit. 

In this Chapter we will critically discuss the following topics 

— Section [2] The realistic evaluation of the total accuracy in the test performed in 
rec ent years with the existing Earth's artificial satellites LAGEOS and LAGEOS 
II (jCiufolini and Pavlislliool ICiufolini et alJbood : iRies et alj|200sh . LAGEOS was 
put into orbit in 1976, followed by its twin LAGEOS II in 1992; they are passive, 
spherical spacecraft entirely covered by retroreflectors which allow for their accurate 
tracking through laser pulses sent from Ear th-based groun d stations according to 
the Satellite Laser Ranging (SLR) technique l|Degnanll 19851 ). They orbit at altitudes 
of about 6000 km (gilageos = 12270 km, «LAGEOS II = 12163 km) in nearly 
circular paths (eLAGEOS = 0.0045, clageos II — 0.014) inclined by 110 deg and 
52.65 deg, respectively, to the Earth's equator. The Lense-Thirring effect for their 
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nodes amounts to about 30 milliarcseconds per year (mas yr _1 ) which correspond 
to about 1.7 m yr -1 in the cross-track directioro at the LAGEOS altitudes. 
The idea of measuring the Lense-Thirring node rate with the just launched LA- 
GEOS s atellite, along with the other SLR targets orbiting at that time, was put 
forth by ICueusi and Proverbiol 1 1978h . Tests have started to be effectively per- 
form ed later by using the LAG EOS and LAGE OS II satellites i Ciufolini et al.l 
1996), according to a strategy by ICiufolinil (|l996T ) involving the use of a suitable 
linear combination of the nodes fi of both satellites and the perigee to of LAGEOS 
II. This was done to reduce the impact of the most relevant source of systematic 
bias, i.e. the mismodelling in the even (I = 2,4,6...) zonal (m = 0) harmonic 
coefficients Jg of the multipolar expansion of the Newtonian part of the terrestrial 
gravitational potential due to the diurnal rotation (they induce secular precessions 
on the node and perigee of a terrestrial satellite much larger than the gravito- 
magnetic ones. The Ji coefficients cannot be theoretically computed but must be 
estimated by suitably processing long data sets from the dedicated satellites like 
CHAMP and GRACE; see Section[2]): the three-el ements combination us ed allowed 
for removing the uncertainties in J2 and J4. In IjCiufolini et al. 1998a ) a » 20% 
test was reported by using th^3 EGM96 Earth gravity model ( Lemoine et al. ill 9981 ); 
subsequent detailed analyses showed that such an evaluation of the total error bud- 
get was overly optimistic in view of the likely unreliable c omputation of the total 
bias due to the even zonals l|lorioll2003l ; iRies et al.ll2003al lbh. An analogous, huge 
underestim ation turned out to hold also for the effect of the non-gravitational per- 
turbations ( Milani et al.|[l987l) like the direct solar radiation pressure, the Earth's 
albedo, various subtle thermal effects depend ing on the physical properties of the 
satellites' surfaces and their rotational state (|Inversi and Vespdll99 j | Vespe|[l999l : 
lLucchesi|[200ll . |2002| . |2003| . 12004 iLucchesi et al.ll2004l ; IRies et al.ll2003ab . which the 
perigees of LAGEOS-like satellites are particular ly sensitive to. As a result, the 
realistic total error budget in the test reported in ((ciufolini et al.|[l998al ) might be 
as large as 60 — 90% or (by considering EGM96 o nly) even more. 
The observable use d bvlCiufolini and Pavlisl (120041) with the GRA CE-only EIGEN- 
GRACE02S model dReigber et al.ll2005t ) and bv lRies et all <|200ct ) with other more 
recent Earth gravity models was the following line ar combination^ of the n odes of 
LAGEOS and LAGEOS II, ex plicitly compute d bv llorio and Moreal (|2004h follow- 
ing the approach put forth bv lCiufolinil 1 19961 ) 



f = n 



LAGEOS 



+ Cl n 



LAGEOS II 



(3) 



3 A perturbing acceleration like Ajjy is customarily projected onto the radial f , transvers e 
t and cross-track O directions of an orthogonal frame comoving with the satellite (Soffcl 1989); 
it turns out that the Lense-Thirring node precession affects th e cross-track component o f the 
orbit according to zIz^lt ~ asiniAQur (eq. (A65), p. 6233 in (Christodoulidis et al. 1988)). 

4 Contrary to the subsequent models based on the dedicated satellites CHAMP 
(http:// www-app2.gfz-potsdam.de/pbl/op/champ/index_CHAMP.html) and GRACE 
(http://www-app2.gfz-potsdam.de/pbl/op/grace/index_GRACE.html), EGM96 relies upon 
multidecadal tracking of SLR data of a constellation of geodetic satellites including LAGEOS 
and LAGEOS II as well; thus the possibility of a sort of a — priori 'imprinting' of the 
Lense-Thirring effect itself, not solved-for in EGM96, cannot be neglected. 

5 See also llPavlis! l2002; Rics et al. 2003___3) . 
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where 



^LAGEOS _ COS tLAGEOS ( 1 ~ e LAGEOS II ^ ( "LAGEOS II \ ^ 



^LAGEOS II C os i L AGEOS II \ 1 - e L,AGEOS / V "LAGEOS 

(6) 

The coefficients Q i of the ahasing classical node precessions 1 Kaula|[l966l ) fi c iass — 
Q f Jp induced b y the even zonals have been analytically worked out up to I — 20 
in, e.g. (|lorioll2003h ; they yield c\ — 0.544. The Lense-Thirring signature of eq. 
amounts to 47.8 mas yr _1 . The combination eq. ((5} allows, by construction, to 
remove the aliasing effects due to the static and time- varying parts of the first even 
zonal J2. The nominal (i.e. computed with the estimated values of Jg, t = 4,6...) 
bias due to the remaining higher degree even zonals would amount to about 10 5 
mas yr _1 ; the need of a careful and reliable modeling of such an important source of 
systematic bias is, thus, quite apparent. Conversely, the nodes of the LAGEOS-type 
spacecraft are directly affected b y the non-grav i tational accelerat i ons at a ~ 1% 
level of the Lense-Thirring effect (|Lucchesill200ll . |2002l . I200A |2004| ; | Lucchesi et al.l 
2004). For a comprehensive, up-to-date overview of the numerous a nd subtle iss ues 
concerning the measurement of the Lense-Thirring effect see, e.g., jlorioll2007bh . 
Section [3] Another approach which could be followed in extracting the Lense- 
Thirring effect from the data of the LAGEOS-type satellites. 

Section[4] The possibility that the LARES mission, recently approved by the Italian 
Space Agency (ASI), will be able to measure the Lense-Thirring node precession 

wit h an accuracy of the order of 1 %. 

In ( Van Patten and Everittlll976al ; rvan Patten and Everitdll976bl ) it was proposed 
to measure the Lense-Thirring precession of the nodes Q of a pair of counter- 
orbiting spacecraft to be launched in terrestrial polar orbits and endowed with 
drag-free apparatus. A somew h at eq uivalent, cheaper version of such an idea was 
put forth in 1986 by ICiufolinil ( 19861 ) who proposed to launch a passive, geodetic 
satellite in an orbit identical to that of LAGEOS apart from the orbital planes 
which should have been displa ced by 18 deg a part. The measurable quantity was, 
in the case of the proposal by ICiufolinil ( 19861 ). the sum of the nodes of LAGEOS 
and of the new spacecraft, later named LAGEOS III, LARES, WEBER-SAT, in 
order to cancel to a high level of accuracy the corrupting effect of the multipoles 
of the Newtonian part of the terrestrial gravitational potential which represent 
the major source o f systematic erro r (see Section Q. Altho ugh extensively studied 
by various groups jRies et al.lll989l; ICiufolini et aL|[l998bh . such an idea was not 
implemented for many years. In (jlorio et alj|2002h it was proposed to include also 
the data fr om LAGEO S II by using a different observable. Such an approach was 
proven m |lorioll2005ah to be potentially useful in making the constraints on the 
orbital configuration of the new SLR satellite less stringent than it was originally 
required in view of the recent improvements in our knowledge of the classical part of 
the terrestrial gravitational potential due to the dedicated CHAMP and, especially, 
GRACE missions. 

Since reaching high altitudes and minimizing the unavoidable orbital injection er- 
rors is expensive, it was explored the possibility of discard ing LAGEOS and LA- 
GEOS II using a low -altitude , nearly polar orbi t for LARES l|Lucchesi and Paolozzil 
l200ll ; ICiufohnj|2006r i. but in (|lorioll2002l . l2007cl ) it was proven that such alternative 
approaches are not feasible. It was also suggested that LARES would be able to 
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probe alternative the ories of gravi t y JCiufolinilbOQ-l ). but also in this case it turned 
out to be impossible (|lorioll2005b1 . l2007dh . 

The stalemate came to an end when ASI recently made the following official an- 
nouncement ( http: / /www. asi.it /SiteEN/MotorSearchFullText.aspx?keyw=LARES ) : 
"On February 8, the ASI board approved funding for the LARES mission, that will 
be launched with VEGAs maiden flight before the end of 2008. LARES is a passive 
satellite with laser mirrors, and will be used to measure the Lense-Thirring effect." 
The italian version of the announcement yields some more information specifying 
that LARES, designed in collaboration with National Institute of Nuclear Physics 
(INFN), is currently under construction by Carlo Gavazzi Space SpA; its Princi- 
pal Investigator (PI) is I. Ciufolini and its scientific goal is to measure at a 1% 
level the Lense-Thirring effect in the gravitational field of the Earth. Concerning 
the orbital configuration o f LARES , In one of the latest communication to INFN, 
Rome, 30 January 2008, l|Ciufolinil l2008al ) writes that LARES will be launched 
with a semimajor axis of approximately 7600 km and an inclination between 60 
and 80 deg. More precise information can be retrieved in Section 5.1, pag 9 of 
the document Educational Payload on the Vega Maiden Flight Call For Cube- 
Sat Proposals, European Space Agency, Issue 111 February 2008, downloadable at 
|http://esamultimedia.esa.int/docs/LEX-EC/CubeSatf o20CFP%20issue%201.pdf. It 
is written there that LARES will be launched into a circular orbit with altitude 
h — 1200 km, corresponding to a semimajor axis aLARES — 7578 km, and inclina- 
tion i = 71 deg to the Earth's equator. Latest informatior0 point towards a launch 
at the end of 2009 with a VEGA rocket in a circular orbit inclined by 71 deg to 
the Earth's equator at an altitude o£| 1450 km corresponding to a semimajor axis 
of a LR — 7828 km. More or less the same has been reported by ICiufolinil l|2008bh 
to the INFN in Villa Mondragone, 3 October 2008. 



2 The systematic error of gravitational origin in the LAGEOS-LAGEOS II 

test 



The realistic evaluation of the total error budget of the LAGEOS-LAGEOS II node test 



Iorioll2005a 



20061 



Ciufolini and Pavlis 2004) raised a lively debate IjCiufolini and Pavlia2005l ; ICiufolini et al 



2007e!;lLucc hcsi 2005), mainly focussed on the impact of the 



static and time- varying parts of the Newtonian component of the Earth's gravitational 
potential through the secular precessions induced on a satellite's node. 

In the real world the path of a probe is not only affected by the relativistic gravit- 
omagentic field but also by a huge number of other competing classical orbital pertur- 
bations of gravitational and non-gravitational origin. The most insidious disturbances 
are those induced by the static part of the Newtonian component of the multipolar 
expansion in spherical har monic jj Jp, £ = 2,4,6,... of the gravitational potential of 



the central rotating mass ( Kaulal 19661 ): they affect the node with effects having the 



same signature of the relativistic signal of interest, i.e. linear trends which are orders of 
magnitude larger and cannot be removed from the time series of data without affecting 



6 See on the WEB http://www.esa.int/esapub/bulletin/bulletinl35/bull35f_bianchi.pdf 

7 I thank Dr. D. Barbagallo (ESRIN) for having kindly provided me with the latest details 
of the orbital configuration of LARES. 

8 The relation among the even zonals Ji and the normalized gravity coefficients CVo which 

-VW+TC m . 



are customarily determined in the Earth's gravity models, is Jt 
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the Lense-Thirring pattern itself as well. The only thing that can be done is to model 
such a corrupting effect as most accurately as possible and assessing the impact of the 
residual mismodelling on the measurement of the frame-dragging effect. The secular 
precessions induced by the even zonals of the geopotential can be written as 



^geopot 



(7) 



where the coefficients Qj,£ = 2, 4, 6, ... depend on the parameters of the Earth (GM 
and the equatorial radius R) and on the semimajor axis a, the inclination i and the 
eccentricity e of the satellite. For example, for I — 2 we have 



3 (R\ 2 _ 

'2 n \a) (1 



(l- e 2) 



1\1 ' 



(8) 



n = y/GM/a 3 is the Keplerian mean motion. They have been analytically computed 
up to I — 20 in, e.g., ( IoriolEooil '). Their mismodelling can be written as 



sn gcopot <^2\n.i\sj e 



(9) 



1=1 



where 5Jg represents our uncertainty in the knowledge of the even zonals Jg 

A common feature of all the competing evaluations so far published is that the 
systematic bias due to the static component of the geopotential was always calculated 
by using the released (more or less accurately calibrated) sigmas crj e of one Earth 
gravity model solution at a time for the uncertainties SJg. Thus, it was said that the 
model X yields a x% error, the model Y yields a y% error, and so on. 

Since a trustable calibration of the formal, statistical uncertainties in the estimated 
zonals of the covariance matrix of a global solution is always a difficult task to be im- 
plemented in a reliable way, a much more realistic and conservative approach consists, 
instead, of taking the differenc^B 



AJt = \Jt(X)-Jt<X)\, 1 = 2,4,6, 



(10) 



of the estimated even zonals for different pairs of Earth gravity field solutions as repre- 
sentative of the real uncertainty S Jg in the zonals l|Lerch et al.lll994h . In Table E-Table 
[9]we present our results for the most recent GRACE-based models released so far by dif- 
ferent institutions and retrievable on the Internet alH http : //icg e m.gfz-potsdam.de/ICGEM/ICGEM.html 
The models used are EIGEN-GRACE02 S (|Reigberet al.ll2005l) from GFZ (Potsdam, 
Germany), GGM02S dTaplev et alj|2005l) and~GGM03S l|Taplev et al.ll2007l ) from CSR 
(Aus tin, Texas), ITG-Grace02s (jMaver-Giirr et allbood ) and ITG-Grace03s (jMaver-Giirrl 
|2007| ) from I GG (Bonn, Germa ny), JEM01-RL03B from JPL (NASA, USA) and AIUB- 
GRACE01S l|jaggi et alj|2008h from AIUB (Switzerland). Note that this approach was 
explicitly followed also by ICiufolinj (|l99rj ) with the JGM3 and GEMT-2 models. In 



-20 



Table[T]-Table|5]we quote both the sum i 
mismodelled terms 

ft = 



of the absolute values of the individual 



ALAGEOS . ALAGEOS II A t 



(11) 



Sec Fig. 5 of l|Luccheaill2007l ) for a comparison of the estimated C40 in different models. 
10 I thank M Watkins (JPL) for having provided me with the even zonals and their sigmas 
of the JEM01-RL03B model. 
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(SAV), and the square root of the sum of their squares yYle=4 ff (RSS); in both 
cases we normalized them to the combined Lense-Thirring total precession of 47.8 mas 

yr- 1 . 



in the eve 
,20, in mas yr 

for the uncertainty in the even zonals we have taken here the difference ACeo 



Table 1 Impact of the 
I^LAGEOS +Cl f)LAGEOS III 



mismodclling 
AJ e , £ = 4,... 



zonal harmonics on fe = 
Recall that J e = -yj2l + 1 C l0 ; 



C 



(X) 



7^(Y) 



between the model X= GGM02S llTaplev et al.l [2005) and the model Y=ITG-Grace02s 
l lMaver-Gurr et al.|[200c^ . GGM02S is based on 363 days of GRACE-only data (GPS and 
intersatellite tracking, neither constraints nor regularization applied) spread between April 4, 
2002 and Dec 31, 2003. The a are formal for both models. ACeo are always larger than the 
linearly added sigmas, apart from I = 12 and I = 18. Values of fe smaller than 0.1 mas yr" 1 
have not been quoted. The Lense-Thirring precession of the combination of eq. J5} amounts 
to 47.8 mas yr - 1 . The percent bias Sfi have been computed by normalizing the linear sum of 
fe , £ = 4, . . . , 20 (SAV) and the square root of the sum of /|, £ = 4, . . . , 20 to the Lense-Thirring 
combined precessions. 



i 


AC eo (GGM02S-ITG-Grace02s) 


0"X + <TY 


h (mas yr 1 ) 


4 


1.9 x 10" 11 


8.7 x 10" 12 


7.2 


6 


2.1 x lO" 11 


4.6 x 10" 12 


4.6 


8 


5.7 x 10" 12 


2.8 x 10~ 12 


0.2 


10 


4.5 x 10~ 12 


2.0 x 10~ 12 




12 


1.5 x 10" 12 


1.8 x 10~ 12 




14 


6.6 x 10~ 12 


1.6 x 10~ 12 




16 


2.9 x 10" 12 


1.6 x 10~ 12 




18 


1.4 x 10" 12 


1.6 x 10~ 12 




20 


2.0 x 10~ 12 


1.6 x 10~ 12 






5fi = 25% (SAV) 


6[i = 18% (RSS) 





Tabl e 2 Bias due to the mismodelling in the even zonals of the models X=ITG-Grace03s 
jMaver-Giirrl 12007ft . based on GRACE-only accumulated normal equations from data out 
of September 2002- April 2007 (neither apriori information nor regularization used), and 
Y=GGM02S llTaplev et al.ll2005l ). The a for both models are formal. AC e0 are always larger 
than the linearly added sigmas, apart from i = 12 and i = 18. 



t 


ACeo (ITG-Grace03s-GGM02S) 


0"X + <TY 


fl (mas yr 1 ) 


4 


2.58 x lO" 11 


8.6 x 10" 12 


9.6 


6 


1.39 x lO" 11 


4.7 x 10" 12 


3.1 


8 


5.6 x lO" 12 


2.9 x lO" 12 


0.2 


10 


1.03 x lO" 11 


2 x 10~ 12 




12 


7 x 10~ 13 


1.8 x 10~ 12 




14 


7.3 x 10~ 12 


1.6 x 10~ 12 




16 


2.6 x 10~ 12 


1.6 x 10~ 12 




18 


8 x 10~ 13 


1.6 x 10" 12 




20 


2.4 x 10~ 12 


1.6 x 10~ 12 






5/i = 27% (SAV) 


<5/i = 21% (RSS) 
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Table 3 Bias due to the mismodelling in the even zonals of the models X = GGM02S 
l lTaplev et al.l l2005h and Y = GGM03S l lTaplev et alj l2007f) retrieved from data spanning 
January 2003 to December 2006. The a for GGM03S are calibrated. AC io are larger than the 
linearly added sigmas for I = 4,6. (The other zonals are of no concern) 



I 


AC m (GGM02S-GGM03S) 




fe (mas yr 1 ) 


4 


1.87 x 10- 11 


1.25 x 10" 11 


6.9 


6 


1.96 x lO" 11 


6.7 x 10~ 12 


4.2 


8 


3.8 x 10~ 12 


4.3 x 10" 12 


0.1 


10 


8.9 x 10~ 12 


2.8 x 10" 12 


0.1 


12 


6 x 10 -13 


2.4 x 10 -12 




14 


6.6 x 10" 12 


2.1 x 10~ 12 




16 


2.1 x 10~ 12 


2.0 x 10~ 12 




18 


1.8 x 10~ 12 


2.0 x 10~ 12 




20 


2.2 x 10" 12 


1.9 x 10~ 12 






Sfi = 24% (SAV) 


Sn = 17% (RSS) 





Table 4 Bi as due to the mism odelling in the even zonals of the mod els X = EIGEN- 
GRACE02S l lReigber et alj|2005f) and Y = GGM03S llTaplev et al.ll2007H . The a for both 
models are calibrated. ACiq are always larger than the linearly added sigmas apart from 
I = 14, 18. 



i 


AC i(i (EIGEN-GRACE02S-GGM03S) 


CT X + CT Y 


ft (mas yr x ) 


4 


2.00 x lO" 11 


8.1 x 10~ 12 


7.4 


6 


2.92 x 10" 11 


4.3 x 10~ 12 


6.3 


8 


1.05 x lO" 11 


3.0 x lO" 12 


0.1 


10 


7.8 x 10~ 12 


2.9 x 10" 12 


0.1 


12 


3.9 x 10~ 12 


1.8 x 10~ 12 




14 


5 x 10~ 13 


1.7 x 10- 12 




16 


1.7 x 10~ 12 


1.4 x 10" 12 




18 


2 x 10~ 13 


1.4 x 10" 12 




20 


2.5 x 10~ 12 


1.4 x 10~ 12 






5/i = 30% (SAV) 


= 20% (RSS) 





The systematic bias evaluated with a more realistic approach is about 3 to 4 times 
larger than one can obtain by only using this or that pa rticular model. The scatte r 
is still quite large and far from the 5 — 10% claimed in ( Ciufolini and Pavlisl booil ). 
In particular, it appears that J4, Jq, and to a lesser extent Jg, which are just the 
most relevant zonals for us because of their impact on the combination of eq. (0, are 
the most uncertain ones, with discrepancies AJg between different models, in general, 
larger than the sum of their sigmas aj e , calibrated or not. 

Another way to evaluate the uncertainty in the LAGEOS-LAGEOS II node test 
may consist of computing the nominal values of the total combined precessions for 
different models and comparing them, i.e. by taking 



E 

1=4 



^LAGEOS _|_ Ci ^LAGEOS II 



[J t (X)-3 t (Y)] 



(12) 



The results are shown in Table [lOl 
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Table 5 Bias due to the mismodelling in the even zonals of t he models X = JE M01-RL03B, 
based on 49 months of GRACE-only data, and Y = GGM03S llTaplev et 3,1.11200711 . The a for 
GGM03S are calibrated. ACiq are always larger than the linearly added sigmas apart from 



£ = 16. 


£ 


AC to (JEM01-RL03B-GGM03S) 


C"X + cry 


fi (mas yr 1 ) 


4 


1.97 x 10" 11 


4.3 x 10~ 12 


7.3 


6 


2.7 x 10~ 12 


2.3 x 10~ 12 


0.6 


8 


1.7 x 10~ 12 


1.6 x 10~ 12 




10 


2.3 x 10" 12 


8 x lO" 13 




12 


7 x 10" 13 


7 x 10~ 13 




14 


1.0 x 10~ 12 


6 x 10~ 13 




16 


2 x 10~ 13 


5 x 10~ 13 




18 


7 x 10" 13 


5 x 10~ 13 




20 


5 x 10~ 13 


4 x 10~ 13 






Sfi = 17% (SAV) 


Sfi = 15% (RSS) 





Table 6 Bias due to the mismodelling in t he even zonals of the models X = JEM01-RL03B 
and Y = ITG-Grace03s llMaver-Gurrll2007ft . The a for ITG-Grace03s are formal. AC m are 
always larger than the linearly added sigmas. 



1 


AC m (JEM01-RL03B-ITG-Grace03s) 


o"X + &Y 


ft (mas yr 1 ) 


4 


2.68 x 10" 11 


4 x 10" 13 


9.9 


6 


3.0 x 10~ 12 


2 x 10~ 13 


0.6 


8 


3.4 x 10~ 12 


1 x 10~ 13 


0.1 


10 


3.6 x 10~ 12 


1 x 10" 13 




12 


6 x 10" 13 


9 x 10" 14 




14 


1.7 x 10" 12 


9 x 10" 14 




16 


4 x 10~ 13 


8 x 10" 14 




18 


4 x 10~ 13 


8 x 10" 14 




20 


7 x 10~ 13 


8 x 10" 14 






<5m = 22% (SAV) 


Sfi = 10% (RSS) 





Table 7 Aliasing effect of the mism odelling in the even zonal harmo nics estimated in the 
X=ITG-Grace03s jMaver-Giirri 12007ft and the Y=EIGEN-GRACE02S jReigber et al.ll2005f) 
models. The covariance matrix a for ITG-Grace03s are formal, while the ones of EIGEN- 
GRACE02S are calibrated. ACiq are larger than the linearly added sigmas for £ = 4, ...,20, 
apart from £ = 18. 



t 


AC eo (ITG-Grace03s-EIGEN-GRACE02S) 


<J X + cy 


ft (mas yr x ) 


4 


2.72 x lO" 11 


3.9 x 10~ 12 


10.1 


6 


2.35 x 10- 11 


2.0 x 10~ 12 


5.1 


8 


1.23 x lO" 11 


1.5 x 10" 12 


0,1 


10 


9.2 x 10" 12 


2.1 x 10~ 12 


0.1 


12 


4.1 x 10~ 12 


1.2 x 10" 12 




14 


5.8 x 10" 12 


1.2 x 10" 12 




16 


3.4 x 10~ 12 


9 x 10" 13 




18 


5 x 10" 13 


1.0 x 10" 12 




20 


1.8 x 10~ 12 


1.1 x 10" 12 






<5 M = 37% (SAV) 


Sfi = 24% (RSS) 
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Table 8 Bias due to the mismodclling in the even zonals of the mode ls X = JEM01-R L03B, 
based on 49 months of GRACE-only data, and Y = AIUB-GRACE01S llJaggi et al.l200Sfl . The 
latter one was obtained from GPS satellite-to-satellite tracking data and K-band range-rate 
data out of the period January 2003 to December 2003 using the Celestial Mechanics Approach. 
No accelerometer data, no de-aliasing products, and no regularisation was applied. The a for 
AIUB-GRACE01S are formal. AC^o are always larger than the linearly added sigmas. 



t 


AC to (JEM01-RL03B— AIUB-GRACE01S) 


<J X + <?Y 


ft ( mas y r 1 ) 


4 


2.95 x 10" 11 


2.1 x 10" 12 


11 


6 


3.5 x 10~ 12 


1.3 x 10~ 12 


0.8 


8 


2.14 x 10~ n 


5 x 10" 13 


0.7 


10 


4.8 x 10~ 12 


5 x 10~ 13 




12 


4.2 x 10 -12 


5 x 10 -13 




11 


3.6 x 10" 12 


5 x 10~ 13 




16 


8 x 10~ 13 


5 x 10~ 13 




18 


7 x 10~ 13 


5 x 10~ 13 




20 


1.0 x 10" 12 


5 x 10" 13 






S/j, = 26% (SAV) 


Sfi = 23% (RSS) 





Table 9 Bi as due to the mism odelling in the even zonals of the models X = EIGEN- 
GRACE02S l lReigber et alj|2005ft and Y = AIUB-GRACE01S llJaggi et al.ll2008f) . The a for 
AIUB-GRACE01S are formal, while those of EIGEN-GRACE02S are calibrated. ACeo are 



larger 


than the linearly added sigmas for I = 4,6, 8, 16. 






i 


ACeo (EIGEN-GRACE02S— AIUB-GRACE01S) 


CTK + CT Y 


ft (mas yr 1 ) 


4 


2.98 x lO" 11 


6.0 x 10" 12 


11.1 


6 


2.29 x lO" 11 


3.3 x 10" 12 


5.0 


8 


1.26 x lO" 11 


1.9 x 10" 12 


0.4 


10 


6 x 10~ 13 


2.5 x 10~ 12 




12 


5 x 10~ 13 


1.6 x 10~ 12 




14 


5 x 10~ 13 


1.6 x 10~ 12 




16 


2.9 x 10~ 12 


1.4 x 10" 12 




18 


6 x 10~ 13 


1.4 x 10~ 12 




20 


2 x 10~ 13 


1.5 x 10~ 12 






Sfi = 34% (SAV) 


Sfj, = 25% (RSS) 





A different approach that could be followed to take into account the scatter among 
the various solutions consists in computing mean and standard deviation of the entire 
set of values of the even zonals for the models considered so far, degree by degree, and 
taking the standard deviations as repr esentative of the u ncertainties SJg,£ = 4, 6, 8, .... 
It yields Sfi — 15%, in agreement with lfUes et alj ( 20081 ) . 

It must be recalled that also the further bias due to the cross-c oupling between Ji 
and the orbit inclination, evaluated to be about 9% in (|lorioll2007el h must be added. 



3 A new approach to extract the Lense-Thirring signature from the data 

The technique adopted so far in (|CiufoUni and Pavlislliooi ) and ijRies et al.ll2008l ) to 
extract the g ravitomagentic signal from the LAGEOS and LAGEOS II data is described 
in detail in |Lucchesi and BalminolkoOfj ; iLucchesjlwrh . In both the approaches the 
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Table 10 Systematic uncertainty 8fi in the LAGEOS-LAGEOS II test evaluated by taking 
the absolute value of the difference between the nominal values of the total combined node 
precessions due to the even zonals for different models X and Y, i.e. fis eo P ot (X) — i?s eo P ot (Y) . 



Models compared 




AIUB-GRACE01S— JEM01-RL03B 


20% 


AIUB-GRACE01S— GGM02S 


27% 


AIUB-GRACE01S— GGM03S 


3% 


AIUB-GRACE01S-ITG-Grace02 


2% 


AIUB-GRACEO 1 S- ITG-Grace03 


0.1% 


AIUB-GRACE01S— EIGEN-GRACE02S 


33% 


JEM01-RL03B— GGM02S 


7% 


ip,mdi RT.n^R— nciMrns 


1 7% 
J- 1 /o 


JEM01-RL03B-ITG-Grace02 


18% 


JEM01-RL03B-ITG-Grace03s 


20% 


JEM01-RL03B— EIGEN-GRACE02S 


13% 


GGM02S— GGM03S 


24% 


GGM02S- ITG-Grace02 


25% 


GGM02S— ITG-Grace03s 


27% 


GGM02S— EIGEN-GRACE02S 


6% 


GGM03S-ITG-Grace02 


1% 


GGM03S- ITG-Grace03s 


3% 


GGM03S— EIGEN-GRACE02S 


30% 


ITG-Grace02-ITG-Grace03s 


2% 


ITG-Grace02— EIGEN-GRACE02S 


31% 


ITG-Grace03s— EIGEN-GRACE02S 


33% 



Lense-Thirring force is not included in the dynamical force models used to fit the 
satellites' data. In the data reduction process no dedicated gravitomagnetic parameter 
is estimated, contrary to, e.g., station coordinates, state vector, satellites' drag and 
radiation coefficients Cjj and Cr, respectively, etc.; its effect is retrieved with a sort 
of post-post-fit analysis in which the time series of the computeof^l "residuals" of the 
nodes with the difference between the orbital elements of consecutive arcs, combined 
with eq. ([5J, is fitted with a straight line. 

In order to enforce the reliability of the ongoing test it would be desirable to proceed 
following other approaches as well. For instance, the gravitomagnetic force could be 
explicitly modelled in terms of a dedicated solve-for parameter (not necessarily the 
usual PPN 7 one) to be estimated in the least-square sense along with all the other 
parameters usually determined in fitting the LAGEOS-type satellites data, and the 
resulting correlations among them could be inspected. Moreover, one could also look 
at the changes in the values of the complete set of the estimated parameters with and 
without the Lense-Thirring effect modelled. 

A first, tentative step towards the implementation of a similar strategy with the LA- 
GEOS satellites in terms of the PPN parameter has been recently taken by Combrinck 
( 20081 ). 



The expression "residuals of the nodes" is used, strictly speaking, in an improper sense 
because the Keplerian orbital elements are not directly measured. 
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4 On the LARES mission 



The combination that should be used for measur ing the Lense-Thirring effect with 
LAGEOS, LAGEOS II and LARES is (|lorioll2005ah 



^LAGEOS + ^LAGEOS II + fc^ARES (13) 



The coefficients fci and k^ entering eq. (|13[) are defined as 

12 LARES ^LAGEOS _ ^LAGEOS ^LARES 
^1 ^ LAGEOS II ^LARES _q^ AUES ^LAGEOS II 



0.3586, 



ALACEOS AL AG EOS 1I_ AL AG EOS II ALAGEOS 

h. ".2 ".4 ".2 J/ .4 fl 07^1 

™* ^LAGEOS II |jLARES_^LARES ^LAGEOS II U.UJtJL. 



(14) 



The combination eq. (|13p cancels out, by construction, the impact of the first two even 
zonals; we have used iilr = 7828 km, ilr = 71.5 deg. The total Lense-Thirring effect, 
according to eq. (fl3|) and eq. (fl4)) . amounts to 50.8 mas yr _1 . 



4.1 A conservative evaluation of the impact of the geopotential on the LARES mission 

The systematic error due to the uncancelled even zonals Jq, Jg, ... can be conservatively 
evaluated as 



ALAGEOS , , ALAGEOS II , , ALARES 



e=6 



SJg (15) 



Of crucial importance is how to assess 5J£. By proceeding as in Section [2] and by 
using the same models up to degree £ — 60 because of the lower altitude of LARES 
with respect to LAGEOS and LAGEOS II which brings into play more even zonals, we 
have the results presented in Table [TT] They have been obtained with the standard and 
widely used Kaula approach (Kaula 1966) in the following way. We, first, calibrated our 
numerical calculation with the analy tical ones p erformed with the explicit expressions 
for tl £ worked out up to I — 20 in (|lorioll2003l ) ; then, after having obtained identical 
results, we confidently extended our numerical calculation to higher degrees by means 
of two different softwares. 

It must be stressed that they may be still optimistic: indeed, computations for 
£ > 60 become unreliable because of numerical instability of the results. 

In Table [12] we repeat the calculation by u sing fo r 8 J/ t he covariance matrix sigmas 
aj t ; also in this case we use the approach by iKaulal (119661 ) up to degree £ — 60. 

If, instead, one assumes 5j£ = s$_, £ — 2,4, 6, ... i.e., the standard deviations of the 
sets of all the best estimates of j£ for the models considered here the systematic bias, 
up to £ = 60, amounts to 12% (SAV) and 6% (RSS). Again, also this result may turn 
out to be optimistic for the same reasons as before. 

It must be pointed out that the evaluations presented here rely upon calculat ions of 
the c oefficients Q g performed with the well known standard approach by Kaula (|Kaulal 
Il966h : it would be important to try to follow also different computational strategies in 
order to test them. 
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Table 11 Systematic percent uncertainty 8fi in the combined Lense-Thirring effect with 
LAGEOS, LAGEOS II and LARES according to cq. JT5} and &J e = AJ e u p to degree 1 = 60 
for the global Earth's gravity solutions considered here; the approach by (Kaula 1966) has 
been followed. For LARES we adopted a LR = 7828 km, i LR = 71.5 deg, e L R = 0.0. 



Models compared (SJ e = AJ e ) Sfi (SAV) Sfj, (RSS) 



AIUB-GRACE01S— JEM01-RL03B 


23% 


16% 


AIUB-GRACE01S-GGM02S 


16% 


8% 


AIUB-GRACE01S-GGM03S 


22% 


13% 


AIUB-GRACE01S-ITG-Grace02 


24% 


15% 


AIUB-GRACE01S-ITG-Grace03 


22% 


14% 


AIUB-GRACE01S— EIGEN-GRACE02S 


14% 


7% 


JEM01-RL03B— GGM02S 


14% 


9% 


JEM01-RL03B— GGM03S 


5% 


3% 


JEM01-RL03B-ITG-Grace02 


4% 


2% 


JEM01-RL03B-ITG-Grace03s 


5% 


2% 


JEM01-RL03B— EIGEN-GRACE02S 


26% 


15% 


GGM02S— GGM03S 


13% 


7% 


GGM02S- ITG-Grace02 


16% 


8% 


GGM02S- ITG-Grace03s 


14% 


7% 


GGM02S— EIGEN-GRACE02S 


14% 


7% 


GGM03S-ITG-Grace02 


3% 


2% 


GGM03S-ITG-Grace03s 


2% 


0.5% 


GGM03S— EIGEN-GRACE02S 


24% 


13% 


ITG-Grace02-ITG-Grace03s 


3% 


2% 


ITG-Grace02— EIGEN-GRACE02S 


25% 


14% 


ITG-Grace03s— EIGEN-GRACE02S 


24% 


13% 



Table 12 Systematic percent uncertainty 8fi in the combined Lense-Thirring effect with 
LAGEOS, LAGEOS II and LARES according to eq. (15) and SJ e = aj e up to degree 1 = 60 
for the global Earth's gravity solutions considered here; the approach by (Kaula 1966) has 
been followed. For LARES we adopted olr = 7828 km, ilr = 71.5 deg, eLR = 0.0. 



Model (SJe = a e ) 


Sfj, (SAV) 8n (RSS) 


AIUB-GRACE01S (formal) 


11% 


9% 


JEM01-RL03B 


1% 


0.9% 


GGM03S (calibrated) 


5% 


4% 


GGM02S (formal) 


20% 


15% 


ITG-Grace03s (formal) 


0.3% 


0.2% 


ITG-Grace02s (formal) 


0.4% 


0.2% 


EIGEN-GRACE02S (calibrated) 


21% 


17% 



4.2 The impact of some non-gravitational perturbations 

It is worthwhile noting that also the impact of the subtle non-gravitational pertur- 
bations will be different with respect to the original proposal because LARES will 
fly in a lower orbit and its thermal behavior will probably be different with respect 
to LAGEOS and LAGEOS II. The reduction of the impact of the thermal accelera- 
tions, like the Yarkovsky-Schach e ffects, should h ave been reached with two concentric 
spheres. However, as explained by lAndrej (|2007h . this solution will increase the float- 



ing potential of LARES because of the much higher electrical resistivity and, thus, the 
perturbative effects produced by the charged particle drag. Moreover, the atmospheric 
drag will increase also because of the lower orbit of the satellite, both in its neutral 
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and charged components. Indeed, although it does not affect directly t he node fl, it 



induc es a secular decrease of the inclination i of a LAGEOS-like satellite ijMilani et al.l 
Il987h which translates into a further bias for the node itself according to 

in which Si accounts not only for the measurement errors in the inclina tion, but als o 
for any unmodelled/mismodelled dynamical effect on it. According to ( Ioriol [2008bl ) . 
the secular decrease for LARES would amount to 

(!} LR ^°- 6masyr " (17) 

yielding a systematic uncertainty in the Lense-Thirring signal of eq. (|13[1 of about 
3 — 9% yr" 1 . An analogous indirect node effect via the inclin ation could be induced by 
the thermal Yarkovski-Rubincam force as well (|lorioll2008bh . Also the Earth's albedo, 



with its anisotropic components, may have a non-negligible effect. 

Let us point out the following issue as well. At present, it is not yet clear how the 
data of LAGEOS, LAGEOS II and LARES will be finally used by the proponent team 
in order to try to detect the Lense-Thirring effect. This could turn out to be a non- 
trivial matter because of the non-gravitational perturbations. Indeed, if, for instance, 
a combination^ 

^LARES + H ^LAGEOS + ^ ^LAGEOS II ( lg) 

was adopted instead of that of eq. (|13[) , the coefficients of the nodes of LAGEOS and 
LAGEOS II, in view of the lower altitude of LARES, would be 

ALAGEOS II o LARES _ OLARES nLAGEOS II 

/ -- -- ■-- -j 3215 



1 ^lares^lageos n_f2 LAGEOS iif2 LAGEOS 

(19) 

ALARES ALAGEOS nLAGEOS ALARES 

h„ — 5£a iLi ~ J /-a Hi, — 4 7744 

f2 LAGEOS fi LAGEOS n — 17 LAGEOS n 17 LAGEOS •*•'<**• 

and the combined Lense-Thirring signal would amount to 676.8 mas yr" 1 . As a conse- 
quence, the direct and indirect effects of the non-gravitationaP^l perturbations on the 
nodes of LAGEOS and LAGEOS II would be enhanced by such larger coefficients and 
this may yield a degradation of the total obtainable accuracy. 



5 Conclusions 

The so far published evaluations of the total systematic uncertainty induced by the 
even zonal harmonics of the geopotential in the Lense-Thirring test with the combined 
nodes of the SLR LAGEOS and LAGEOS II satellites are likely optimistic. Indeed, 
they are all based on the use of the covariance sigmas, more or less reliably calibrated, 
of the covariance matrices of various Earth gravity model solutions used one at a time 
separately in such a way that the model X yields an error of x%, the model Y yields an 

12 The impact of the geopotential is, by construction, unaffected with respect to the combi- 
nation of eq. (1131) , 

13 The same may hold also for time-dependent gravitational perturbations affecting the nodes 
of LAGEOS and LAGEOS II, like the tides. 
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error y%, etc. Instead, comparing the estimated values of the even zonals for different 
pairs of models allows for a more conservative evaluation of the real uncertainties in 
our knowledge of the static part of the geopotential. As a consequence, the uncertainty 
in the Lense-Thirring signal is about 3 — 4 times larger than the figures so far claimed 
(5-10%), amounting to various tens percent (37% for the pair EIGEN-GRACE02S and 
ITG-GRACE03s, about 25 - 30% for the other most recent GRACE-based solutions). 

Concerning the extraction of the Lense-Thirring signal from the data of the LAGEOS- 
type satellites, different approaches with respect to the one followed so far should be 
implemented in order to do something really new. For instance, the gravitomagnetic 
force should be explicitly included in the dynamical force models of the LAGEOS satel- 
lites and an ad-hoc parameter should be estimated in the least-square sense in addition 
to those determined so far without modelling the Lense-Thirring effect. Moreover, also 
the variation of the values of all the other estimated parameters with and without the 
gravitomagnetic force modelled should be inspected along with their mutual correla- 
tions. 

Applying the strategy of the difference of the estimated even zonals to the ongoing 
LARES mission shows that reaching a 1% measurement of the Lense-Thirring effect 
with LAGEOS, LAGEOS II and LARES maybe difficult. Indeed, since LARES will 
orbit at a lower altitude with respect to the LAGEOS satellites, more even zonal 
harmonics are to be taken into account. Assessing realistically their impact is neither 
easy nor unambiguous. Straightforward calculations up to degree t — 60 with the 
standard Kaula's approach yield errors as large as some tens percent; the same holds 
if the sigmas of the covariance matrices of several global Earth's gravity models are 
used. Such an important point certainly deserves great attention. Another issue which 
may potentially degrade the expected accuracy is the impact of some non-gravitational 
perturbations which would have a non-negligible effect on LARES because of its lower 
orbit. In particular, the secular decrease of the inclination of the new spacecraft due 
to the neutral and charged atmospheric drag induces an indirect bias in the node 
precessions by the even zonals which, in the case of LARES, should be of the order of 
«3-9%yr- 1 . 
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